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^ ' Abstract. We give uniqueness theorem and reconstruction algorithm for the non- 

linearized problem of finding the dielectric anisotropy / of the medium from non-overdeter- 
mined polarization tomography data. We assume that the medium has uniform background 
parameters and that the anisotropic (dielectric permeability) perturbation is described by 
symmetric and sufficiently small matrix-function /. On a pure mathematical level this 
' (-('■ article contributes to the theory of non-abelian Radon transforms and to iterative methods 
Ph' of inverse scattering. 

^ \ 1. Introduction 

We consider the system 



< 



o 

0^ 



9d,T] = TVef{x)rj, x G M^ 6* G §^ (1.1) 



where 



> 
(N 



^^- = E^^a^' (1-2) 



r] at fixed ^ is a function on M with values in Zg, (1-3) 



O ■ Ze = {z eC^ : ze = 0}, (1.4) 

/ is a sufficiently regular function on M"^ with values in A^3 3 

rN ' / ' (1-5) 

^ I (that is in 3x3 complex matrices) with sufficient decay at infinity, 

1T0 is the orthogonal projector on Zq. (1-6) 

In (1.1) the unit vector 9 is considered as a spectral parameter. 

System (1.1) arises in the electromagnetic polarization tomography and is a system 
of differential equations for the polarization vector-function 77 in a medium with zero con- 
ductivity, unit magnetic permeability and appropriately small anisotropic perturbation of 
some uniform dielectric permeability. This anisotropic perturbation of the dielectric perme- 
ability tensor corresponds to the matrix- function / of (1.1). In addition, by some physical 
arguments, / must be skew-Hermition, fij = —fji. For more information on physics of 
the electromagnetic polarization tomography see [Shi], [NS], [Sh3] and references therein 
(and, in particular, [KO] and [A]). 
Let 

cc;G§^ 6'G§i, e^=uxe, (1.7) 
Hi = -quj, 112 = for 7] G Zq, (1.8) 
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where 

§1 = {61 G : eio = 0}, (1.9) 

X denotes vector product, Ze is defined by (1.4). 
From (1.1)-(1.9) it follows that 

where /x is related with r] of (1.1) by (1.8) and is a C^-valued function on M.^ for fixed uj 
and 6. 

We consider also (1.10) for fj, taking its values in Al2,2 (that is in 2 x 2 complex 
matrices). Let denote the solution of (1.10) such that 



lim n+{x + se,e,io) ^ Id for x eW, (1.11) 

s— »— oo 

where Id is the 2x2 identity matrix. Let 

s{x,e,uj)^ lim n+{x + se,e,uj), {x,e)eT^, (1.12) 

s— »+oo 



where 

T^ = {{x,e): xeXe, eeSl}, a; g S^, 
X0 = ReZ0, ^GS^ 



(1.13) 



where is defined by (1.9), is defined by (1.4). In addition, /i+ and S are well defined 
due to (1.5). 
Note that 

c T§^ u G §^ (1.14) 

where 

TS"^-! = {(x, 6*) : a; G M'^, ^ G S^-\ = 0}. (1.15) 

In addition, we interpret TEi^~^ as the set of all rays in R'^. As a ray 7 we understand a 
straight line with fixed orientation. If 7 = {x,e) G TE>'^~\ then 

7 = {y G : y = X + s9, s G M} (modulo orientation) and 9 gives the orientation of 7. 
Note also that 

r^RiR^x§\ (1.16a) 

or, more precisely, 

{x,9)eT^'^x = Cid^ + C2UJ, e = (6,6)eK', 9eSl^S\ (1.166) 
where a; G S^, = a; x ^. In addition, we consider (^, ^) G x as coordinates on T^. 
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One can see that S of (1.12) is a matrix-function on 



{(7,0;): 7 = (x,^)eT§2, ooeSl}, 



(1.17) 



where E>1 is defined as in (1.9). On the other hand, one can show that 

S{x,9,uj) at fixed a; G and {x,9)eTS^ 

uniquely determines S{x,9,-) on §^ and, (1-18) 
as a corollary, S can be considered as a matrix — function on T§^. 

The matrix-function S can be considered as a non-abelian ray transform of /. See 
[MZ], [V], [Sh2], [N], [FU], [E], [M], [DP], [P] and references therein for some other non- 
abelian ray transforms. 

In the present work we say that S of (1.11)-(1.13) is the polarization ray transform of 

/• 

Using the terminology of the scattering theory one can say also that S is the "scat- 
tering" matrix for system (1.10). 

The basic problem of the polarization tomography in the framework of the model 
described by (1.1), (1.10) consists in finding / on from S on A, where A is some 
appropriate subset of S of (1.17). It is especially natural to consider this problem for the 
case when dim A = 3. 

From results of [NS] it follows that there is a non-uniqueness in this problem if / is 
not symmetric even if S is given on A = E. Results of [NS] also imply a local uniqueness 
theorem (up to a natural obstruction if / is not symmetric) for the case when S is given 
on S (or on in the sense (1.18)). 

In the present work we consider the following inverse problem for equations (1.1), 
(1.10). 

Problem 1.1. Find symmetric /, = fji, from S on A (or from partial information 
about S on A), where 

A={(7,u;): 7er^, 00 E {u;\ . . . ,00''}}, (1.19) 

where co'^ , . . . , co'^ are some fixed points of 

One can see that Problem 1.1 is a version of the aforementioned basic problem of the 
polarization tomography with dim A = 3, see definitions (1.13), (1.17), (1.19). 

The main results of the present work consist in uniqueness theorem and reconstruction 
algorithm for nonlinearized Problem 1.1 with sufficiently small /, where only the element 
Sii of S = (Sij) on A is used as the data and where 

k = 6, uj^ = ei, a;^ = 62, UJ^ = 63, ^ 
a;4 = (ei + e2)/\/2, co^ = (d + es)/V2, a;^ = (e2 + 63)7^2, 

where ei, 62, 63 is the basis in R^. See Sections 2, 3 and, in particular. Theorem 3.1. 
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One can see that this reconstruction is non-overdetermincd: wc reconstruct 6 functions 
fiji 1 < i < J < 3, on from 6 functions u), u G {w^, . . . , uj^}, of 3 variables. 

Our reconstruction is iterative and its first apprximation more or less coincides with 
the linearized polarization tomography reconstruction of Section 5.1 of [Shi]. In addition, 
we give estimates on the reconstruction error / — /" for the approximation /" with number 
n e N, see Theorem 3.1 . To our knowledge even f — was not estimated rigorously in 
the literature. 

The main results of the present work are presented in detail in Sections 2 and 3. 
Some possible development of the present work and some open questions are mentioned 
in Section 6. 

2. Reconstruction algorithm 

Consider the classical ray transform / defined by the formula 

if{7) = J f{x + se)ds, j = {x,e)eTS''-\ (2.1) 



for any complex-valued sufficiently regular function / on with sufficient decay at infinity, 
where TE>^~^ is defined by (1.15) (and where d — 2 oi d = 3). 

We use the following Radon- type inversion formula for / in dimension d = 2: 

/(^) = ^/ h'{x0^,0)d0, h'{s,0) = j-h{s,0), 

(2.2) 

h{s,0)^-p.v. I ^-^dt, 

TT . S-t 



where g{s,0) = If(s0^,0), x = ixi,X2) G R^ ^ = (01,02) e §\ 0^ = i-02,0i), seR,d0 
is the standard element of arc length on S"*^. 

We use the following slice by slice reconstruction of / on from g = If on T^^ of 
(1.13) for fixed to G S^: 

(2.2) 



/ 

tS\y) 



(2.3a) 

Y 
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for each two-dimensional plane Y of the form 

Y = X{^l)+y, |/eX^(§i), (2.36) 

where §^ is defined by (1.9), ^(S^) is the linear span of si, in R?, X^(§^) is the orthogonal 
complement of -^(S^) in M^, TS"'^(y) is the set of all oriented straight lines lying in Y . In 
addition, 

= u,ex-(§i)^^'(^(si) + v)- (2-4) 
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fjj 




/l2 


= 5'u;4 


/l3 




/23 





Consider the three-dimensional transverse ray transformation J defined by the formula 
(see Section 5.1 of [Shi]): 

J ojf{x + s9)iods, (7,0;) e S, 7 = {x, 
R 

for any TWa^a-valued sufficiently regular function / on R^ with sufficient decay at infinity, 
where ujfuj is defined as in (1.10), E is defined by (1.17). 

We use the following reconstruction of symmetric / on R^ from J/ on A of (1.19) for 
uj^,. . .,u;'^ given by (1.20): 

J = 1,2,3, 

- 2^h^ + /22), 

1 (2.6) 

- 2^fl-l + /33), 



where g^o = Jf\-p and ^ denotes the slice by slice reconstruction via inversion formulas 
(2.2), (2.3) for / from data on F^. 

Now we are ready to present our iterative reconstruction of sufficiently small, sym- 
metric and compactly supported / from the element of S = (Sij) on A, where S is 
defined by (1.12), A is defined by (1.19), (1.20). 

Thus, in addition to (1.5), we assume that 

/ is symmetric, fij = fji, 

f{x) = for \x\ > ro, (2.7) 
/ is sufficiently small. 

Let 

A° = (5n-l)|^, (2.8) 
f = xJa 'AO, (2.9) 

where J^^ denotes the reconstruction via inversion formulas (2.6) for J from data on A, 
X denotes the multiplication by smooth x such that 

xix) = 1 for \x\ < rn, 

- n f ^ (2.10) 

X{x) = tor > ri, 

where ro is the number of (2.7), ri > tq. 
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In our iterative reconstruction, is the first approximation to /. 
From the approximation with number n the approximation f^~^^ with number 
n + 1 is constructed as foUows: 

(1) We find the element //"j^ of //'^+ = (a*^"^) on 

V = {{x,e,u;) : X eR^, 9 eSl, u; e {uj\...,u;^}}, (2.11) 

where satisfies (1.10), (1.11) with /" in place of / in (1.10); 

(2) We find 

S'^^{x,0,u)= lim fj,1+(x + se,e,u;), (x,e,uj) e A, (2.12) 

s— >+oo 

A- = {S,^ - S^,)\^; (2.13) 

(3) Finally, we find 

r+' = x(r + Ji'A-), (2.14) 

where J^^ and x c^re the same that in (2.9). 

Note that in (2.9), (2.14) we do not assume that A'^, A"^ are in the range of J. 
However, J^^g is well-defined on the basis of (2.6) for any 

^= (^^i,...,5r^6), where 
g^ji is a complex — valued sufficiently regular 
function on T^i with sufficient decay at 
infinity (see (1-16)) for each z e {1, . . . , 6}. 



3. Convergence 

We consider 



L°-^-{W') = {u: ueL^iW"), < +oo}, a > 0, (3.1) 



where 

u 



(p) = (i-)^ J e'P''u{x)dx, p e M^ (3.2) 



Ml^,^(R') = ess sup (1 + \p\y\u{p)\. (^-^^ 

We consider 

C^'^'iR^) = {u: ue C["](m2), ||«||(5„,^(K3^ < +oo}, a>0, a>0, (3.4) 
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where u is defined by (3.2), Ct'*' denotes [aj-times continuously difFerentiable functions, 
[a] is the integer part of a, 

ll"llc«-(M^) = M\c-,^(M.'y (3-5) 

ll^llc«.-(R') = (1 + for (3.6) 

|j|<[a], peR' 

\mr<-,<T(W^^ ='^^iNi,N2) for a > [a]. 



^ |9'^'u(p') — 9'^{i(p)| 



^1 = McM,^{R') 



(3.7) 



|JN[a], peM^ p'eM.\ \p-p'\<i \P P \ 



where 



d-'u{p)= ^ J J={J,,J2,J3)e{NiJ0r, \J\=J, + J2 + J3. (3.8) 



In addition, in (3.1)-(3.8) we assume that u, u are A^ni,n2"Valued functions, in general, 
where Alni,n2 is the space of ni x n2 matrices with complex elements, 

|M| = max \Mij\ for M e Mn,,n2- (3-9) 

l<i<ni 

In addition, in the present work we always have that 1 < ni < 3, 1 < 77,2 < 3. 

Lemma 3.1. Let u e L~''"(K^), v e (^"'^(R^), w;/iere a > 0, o" > 3. In addition, 
in general, we assume that u is M.m,n2~ valued and v is M.rm,m2~ valued, where 777.2 = ni 
or/and 712 = rni {and where 1 < ni, n2, mi,m2 < 3). Let 



(3.10) 



(3.11) 



either w = vu for m2 = ni 
or w = uv for n2 = toi. 

Then for each ofw of (3.10) the following estimate holds: 

w e C"*'''(]R^) 

for some positive Ai = Ai(q!, a). 

Lemma 3.1 is proved in Section 4. 
We consider 

L°°^^A) = {g: geL°°{A), |b|||,oo,.(A) < +00}, a > 0, (3.12) 
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where g is complex-valued, A is defined by (1.19) with a;^, . . .^uj^ given by (1.20), 

9= (fl'a;i,---,fi'a;fe). 9= {9uj^ , ■ ■ ■ , 9w'^) , 

9u>^ — 9\r ■' 9lo^ = 9\-n .) 

^ (3.13) 

9oji{p,0)={—f e'P''g^i{x,e)dx, {p,e) eT^i, i = l,...,k, 
Xe 

II^IIl°°.<'(a) ~ II^IIl°°.°-(a)) 

||^||l°°.<-(A) = max ess sup {I + \p\Y\guji{p,d)\, (3-14) 
ie{i,...,fe} (p,0)er^i 

where and Xq are defined according to (1.13). Actually, in (3.12) we consider L°°(A) 
as 

L°°(A) = L°°(r^i) © ... © L°°(r^O- (3.15) 



We assume that 



/ is a Al3,3 — valued function on M^, 

/eL°°''"(R3) for some cr > 3, 
/ is symmetric, fij = /^j, 
/(x) = for |a;| > tq. 



(3.16) 
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X is a nonnegative real — valued function on 
X e C"'(]R^) for some m e N, m > a, 

x{x) > xiv) if \y\ > (3.17) 

x(a;) = 1 for |x| < tq, 
x(a;) = for \x\ > ri, 

where ro, ri are some fixed real numbers, tq < ri. Properties (3.16), (3.17) imply, in 
particular, that 

Xf = f, 

X e C'"''"(M'') for any a > 0. 



3^ . . . (3-18) 



Theorem 3.1. Let f and x satisfy (3.16), (3.17). Let 

<e< £o(«,cr, p). 



(3.19) 

llx||(5i+«,.(Ma) < P 

for some a e]0, 1[, p > 0, eo > 0, where Eq = eo{oi,a,p) is sufficiently small. Let S be the 
polarization ray transform of f (see Section 1 and, in particular, formulas (1.11), (1.12)). 
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Let A be defined by (1.19), (1.20). Then the element Sn of S = (Sij) on A uniquely 
determines f by the iterative reconstruction algorithm of Section 2; in addition, 

/ in L°°''"(m2) as n^+oo, 
ll/-rillco..(M3^<a(«,^,P)(%,^,P))"-'£'^+\ ^" ' 

for some a and b, where f^, n G N, are defined by (2.9), (2.14). 

Theorem 3.1 follows from Lemma 3.1, properties (3.18), and Propositions 3.1, 3.2, 3.3. 

Proposition 3.1. Let 

g e L°°'^(A) for some a > 3, (3.21) 

where g is complex-valued (and A is defined by (2.19), (2.20)). Let J~^^ be defined as in 
(2.9). Then 

J7^g G L°°''^(M^), 

II -1 (3-22) 

for some positive A2. In addition, if g = Jf, where f satisfies (3.16), then g satisfies (3.21) 
and 

Jl^g = /• (3.23) 
Proposition 3.1 is proved in Section 4. 

Proposition 3.2. Let the assumptions of Theorem 3.1 be fulfilled. Let S^Jf on A 
and 6^f on be defined by the formulas 

A° = Jf + d^Jf on A, (3.24) 
J^^A° = / + 5V on M^ (3.25) 

where A°, J^^A^ are defined as in (2.8), (2.9). Then 



5^Jf G L~'"(A), 
S^f e l°°'<"(r3), 

II^Vllroo.arM^ < A3 (a, cT,p)£^ 



(3.26) 



(3.27) 



for some positive A3 = A3 (a, a, p). 

Proposition 3.2 is proved in Section 4. 

Proposition 3.3. Let 

r = x(/ + <^"/), (3.28) 
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where /, x satisfy (3.16), (3.17), 

< £ < ei{a,a,p), 

(3.29) 



ll/ + <^''/llf,oc..r^3^ < £ < £i(a,(7,p). 



/or some o: e]0, 1[, p > 0, £i > 0, where £\ = Si{a, ex, p) is sufficiently small. Let f^'^^ be 
constructed from Sn {for f) and f^ as described in Section 2. Then 

r+' = x(/ + 5"+V), (3.30) 

where 



f3 31) 

ll^'^+Vllf,o.,.,^3^<A4(«,a,p)£||<5-/||. ,^3^ ^ ■ ^ 



/or some positive A4 = A4(q;, (J, p). 

Proposition 3.3 is proved in Section 5. 

To obtain Theorem 3.1 we assume that eo of (3.19) and Si of (3.29) are so small that 

Sq + < El, A4£i < 1, (3 32) 

A4(l + A2A3£o) < beo < 1, 

for some b = b{a,a,p), where A2, A3, A4, £1 are the constants of Propositions 3.1, 3.2, 
3.3. Under these assumptions. Theorem 3.1 follows directly from Propositions 3.2, 3.3 and 
Lemma 3.1. In addition, we use Proposition 3.3 with e given as e + A2A3£^ in terms of e 
of Theorem 3.1. 

4. Proofs of Lemma 3.1 and Propositions 3.1 and 3.2 

Let 

L~'-(M3) = {^ieL°°(M=^): < +00}, a>0, (4.1) 



where || • ||^oo,o-(-]R3^ is defined as in (3.3), 



C"''^(M3) = |^eC'[«](M3) . ||ii||^^^^^^3^ < +00}, a>0, (j>0, (4.2) 
where C^"^ and || • ||(^a,o-(]R3^ are defined as in (3.4) and (3.6), (3.7), 

L°°''^(A) = {^eL°°(A): ||^|Uo.,.(A) < +00}, a>0, (4.3) 

where || • ||l°°><^(A) is defined as in (3.14). In addition, we assume that u of (4.1), (4.2) is 
matrix- valued (of some fixed size), in general, and g of (4.3) is complex- valued. 
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Proof of Lemma 3.1. We use, in particular, that 

U1U2 =Ul* U2-, 

ui*U2{p) — J ui{p — p')u2{p')dp' , peM.^, i'^-'^) 

where ui, U2 are test functions on M^, « is defined by (3.2). In addition, ^i, U2 are 
matrix- valued, in general, where the matrix product is defined in the standard way. 
By the assumptions of Lemma 3.1 we have that 

u e L°°''"(E^), V e C"''"(M^). (4.5) 

Formulas (3.11) follow from (4.4), (4.5), where we use, in particular, that 

< for a > 3, (4.6) 



(4.8) 



{i + \p- p'\y{i + \p'\y ~ (1 + \p\y 

for some positive ci = ci{a). 
Lemma 3.1 is proved. 

Proof of Proposition 3.1. We use that I^^g^j of (2.3), (2.6) can be defined also as 
lr,^Quj = Un,^ where 

{ATT)-\g^{p,d) + g^{p,-d)) = uM. peXe, ee^l, 

where u,^ and -u^^ are related by (3.2), g^^^ and g,^ are related as in (3.13). Indeed, (4.7) 
means that 

(27r)-3y" e'^''g'J"^{x,d)dx={2'K)-^ j e'^''u^{x)dx, 
p e Xe\{x = suj,se R}, Oe^l, u e §^ 

where 

g!y^ix,e)^^{g^{x,e)+g{x,-e)), {x,e)eT^. (4.9) 
Representing p in (4.8) as 

p = pi9-^+P2LV (4.10a) 
and representing x in the left hand side of (4.8) as 

x = ei^^ + 6^ (4.106) 
and integrating (4.8) with e~*P2^2, one can see that u^, of (4.7) is the function such that 

(27r)-2 J e^^^P^g^y^ii^e^ + ^c^, ^)5(6 - i2)diidi2 = 

(4.11) 

(27r)-2 / e'P^^ ''u^{x)6{ojx - Qdx, pi G K\0, G K. 
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It remains to note that it is actually well-known that the determination of itcu|y from 
9uj\rpgi(^Y) (4-9)) (4-11) is equivalent to such a determination on the basis of (2.3), 
where y = ^2UJ. 

Formulas (4.7) imply that 

e l~''^(r3), 

ll«u;|lioo..(M3) < (47r)-^||^a;||L-.-(r„) (4.12) 

if G L°°''^(r^), 

where L°°'^{T^) is considered as L°°'^{A) (see (4.3)) for the case when A is reduced to the 
single r^^. 

Proposition 3.1 follows from (2.6) and (4.12). 

Proof of Proposition 3.2. We consider I0, Dg defined by 

Ieu{x) = J u{x + s9)ds, x e Xg, 9eS^, (4.13) 
R 

+00 

D0u{x) = J u{x + se)ds, a; e R^ de §^ (4.14) 


where tt is a matrix- valued test function on M^, Xq is defined in (1.13). 
We use that 

l_i+{-,e,uj) = Id + D.e{F{-,e,uj)^+{-,e,uj)) on M^ (4.15) 

S{;e,io) ^ Id + l0{F{;e,io)ii+{;e,io)) on Xe, (4.16) 

where fi+ , F, S are defined in Section 1 (see (1.10), (1.11), (1.12)), 9 e S^, u e §^ In 
addition, (4.15) is an integral equation for (4.16) is a formula for S. 

Lemma 4.1. Let u e C^'^i^^) for some cr > 3. Then 



hueC^^^Xe), 

¥eu\\^o,a^Xe) - 27r||w||(5o.a(]R3): ^ ^ 



(4.17) 



where 



C^'-(Xe) = {^7 : geC{X0), \\g\\co,.^Xe)<+^}^ (4-18) 
9{p) = J e'P-g{x)dx, p e Xg, (4.19) 



\9\\cO'-(Xe) - ll^llco.-(Xe), 

l^llco.-(Xe) = sup (1 + \p\y\g{p)\. ^^"^^^ 
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Lemma 4.1 follows from the formula 



(27r)-i^ = «|^^ for g^hu (4.21) 



and definitfons (3.2)-(3.6), (4.18)-(4.20). 
Lemma 4.2. Let 



v{x,d,Uj) = [ n± \ \ n± V \n± > (4-22) 



where u is a Aiz^z- valued function on R"*, X e uj, 9, 6 are vectors of (1.7). Then 

ll^(-,^,^)ll£oc..(R3) < C2||«||^^,.(R3) for u e L-'-(K^), a > 0, (4.23) 
lb(-,^,^)llc«,.(M3) < C2||ti||^„,.(K3^ /or « e ^"'"(R^), a > 0, a > 0, (4.24) 

where C2 is some positive constant. 

Lemma 4.2 follows from definitions (4.22), (3.1)-(3.8). 

Lemma 4.3. Let u e C'«''"(M^) for some a e]0, 1[, cr > 3. Let v e C'i+"''^(M^). Then 

vD-0U e C'"''^(R^), 

< C3(a,(T)||^;i|^i+„,,^^3^||w||^„^^^^3^, 6'e§^ 

/or some positive C3 = 03(0:, cr). 

In Lemma 4.3, m, v are matrix-valued, in general, where the matrix product is defined 
in the standard way. 

Proof of Lemma 4-3. We use that 

D_0u{x) = Gj * u{x) = J Gj{x- y)u{y)dy, 

R" 

Gj{x) = s{ojx)s{e^x)h{ex), 

where u>, 9, 9-^ are related as in (1.7), 6 is the Dirac function. 



(4.26) 



h{s) = 1 for s > 0, 
h{s) = for s < 0. 



(4.27) 



Further in this proof we assume for simplicity that 

9 = ei, 9^ = €2, a; = 63, (4.28) 

o 

where ei, 62, 63 is the basis in R . 
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We use that 

u{Tu2 = {27r)^uiU2, (4.29) 

where ui, U2 are the same that in (4.4). 

We use also that, under assumptions (4.28), 



G+{p) = (27r)"^y e'P''S{x3)6{x2)h{xi)dx = 

w -11 i (^-'^^ 

(27r)-2— / e'P''''h(xi)dxi = (27r)-2 — — = (27r)-^ , p e 

R 

Due to (4.26)- (4.30), we have that 

^(p) = ^^, peK'- (4.31) 

P\ + ^U 

Due to (4.4), (4.31), we have that 

vDeu{p) = J v{j>-p) TI^^q ^/, peR^- (4.32) 
R' 

To prove Lemma 4.3, it is sufficient to prove that 

C3,i(a,a)||i)||^.,^^3 ||^||^„,^^3 

l^i^e«(p)l < {i + \p\Y ' ^ 

| a^i^n(p) | ^ C3,2(Qi,CT)||i)||^,^,,,^^3^||^||^„,,^^3^ 



for p e R^, J = 1, 2, 3 and some positive cs^i, 03^2- 
Proceeding from (4.32) we have that 



vDeu{p)=(^j + J y(p-p')^^dp' = Aip) + Bip), (4.34) 
b'lKi |p'il>i ^ 

v{p-p')u{p)dp' +p.v. J v{p - p) ' ^; dp' = Ai(p) + A2(p)j(4.35) 



^^i/l + IP-p'D'd+lp-l)' (4.36) 

C3,l,l('^)ll^llco.-(R^)ll^llco.- rR^(l+W)"^ 
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where p G M^, C3,i,i, 03,1,2 are some positive constants, ci is the constant of (4.6). 

Estimate (4!33a) follows from (4.34)-(4.38). Estimate (4.33b) follows from (4.34)- 
(4.38) with v{p) replaced by dv{p)/dpj. 

Lemma 4.3 is proved. 

We continue the proof of Proposition 3.2. 

Using (4.15), the property that 

F{;9,iv) = xF{-,e,u), (4.39) 

where x is the function of (3.16), (3.17), and Lemmas 3.1, 4.3 we obtain that 

+00 



(4.40) 



J-1 



F{; e, a;)/x+(-, e, u) = F(-, e,u) + Y, (4-41) 



where 



Wji;^,'^)\\c^,.(U^) < (c3(Q;,cT)||F(-,^,a;)||^,+„,^^^3^)' '||F(-, ^, a;)||^„,^^^3^ < 



qi{q2y 
where 



^ ' (4.43) 
92 = Xi{l + a,a)c2Cs{a,cr)\\x\\^i+^,^(^3y 



and Ai, C2, C3, e are the numbers of Lemmas 3.1, 4.2, 4.3 and Theorem 3.1. 
Due to (4.41)-(4.43), we have that 

||F(-,0,c^)Ai+(-,^,c^)-F(.,^,cc;)||^„,,^^3^< 



+00 5 



E/ yi(/2t (4.44) 
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under the condition that q2S < 1. 

Due to (4.16), (4.44) and Lemma 4.1, we have that 



||5(-,^,a;)-/d-/,F(-,^,a;)|bo,.(x,) < (4-45) 



and, in particular. 



||^n(-,^,a;)-l- Wa;||^o,.(x,) < (4.46) 

-L 

To obtain (3.26) we use (4.46), (4.43) and the following estimate 

^ C4(Q!,(j)||x|lci+«,a(R3)> (4-47) 

where C4 is an appropriate positive constant. In addition, it is assumed that £0 of 
Theorem 3.1 is so small that 

q2So < 1/2. (4.48) 

As a result, proceeding from (4.46) we obtain (3.24), (3.26). 

Finally, (3.27) follows from (3.24)-(3.26) and Proposition 3.1. 
Proposition 3.2 is proved. 

5. Proof of Proposition 3.3 

We will use that in the construction of /"'"'" from and the steps given by (2.13), 
(2.14) can be rewritten as 

r+'=xJX\S^i-l-T^,), (5.1) 
where T{\ is the element of T"^ = (T?J), where 

T"(a;,^,a;)= lim t/"(a; + s^, ^, w), {x,e,u;) e A, (5.2) 

s— »+oo 



dd^U''(x,e,u}) = F''{x,e,u}){ii''+{x,e,oj) - Id), 
lim U''{x + sd,9,u;) = 0, {x,9,oj)eV, 



(5.3) 



where and /x"^"'" are defined as F and of (1.10), (1-11) with in place of /, V is 
defined by (2.11). 

Indeed, (5.2), (5.3) and the definition of /x""'' imply that 

W'ix^e^uj) = ij''+{x,e,uj) ~ Id - D^eF''{x,6,uj), {x,6,Lj)eV, (5.4) 
T^^{x,9,uj) = S^^{x,9,u;)-1-Ieu>r{x)uj, {x,9,Lo)eA, (5.5) 

where I0 and Dq are defined by (4.13), (4.14). Using (5.5), (2.5) one can see that (5.1) is 
equivalent to 

r+' = x^A'(^ii-^ri + <^r)- (s-e) 
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Using that Jj^^Jf^ = (this is completely similar to (3.23)), one can see that f'^'^^ of 
(5.1), (5.6) coincides with of (2.14). 

Thus, it is sufficient to prove Proposition 3.3, where (2.14) is written as (5.1). 

In this proof, in addition to T" and we consider also T and U (defined in a 
completely similar way with T", U"^): 

T{x,e,co)^ lim U{x + se,e,u), {x,e,u)eA, (5.7) 

s— »+oo 



9da;U{x,9,u) = F{x,9,uj){fj.+{x,e,uj) - Id), 
lim U(x + se,e,u;) = 0, (x,e,u;)eV, 



(5.8) 



In addition, (in a completely similar way with (5.4), (5.5)) we have that 

U{x,9,uj) = ii+{x,d,io)- Id- D_0F{x,d,u;), {x,d,u;)eV, (5.9) 

Tn{x,e,u;) = Sn(x,e,oj)-l-l0Ujr(x)u;, (x,d,u;) E A. (5.10) 

One can see that 



xJX\Jf) + xJK\t,, - T-) ^'^'^ xif + <^"+V), 



(5.11) 



where 

5"+V = J-i(Tn-T-). (5.12) 

Thus, to complete the proof of Proposition 3.3 it is sufficient to obtain an appropriate 
estimate on Tn — T{\ (estimate (5.31) given below). 
Let 

rTn=T{\-Tn, (5.13) 

5>+ (5.15) 

Using the definitions of and T one can see that S'^Tn is the element of = ((5'^Ty ), 
where 

S^T = /,(F"(., e, oj)ifi^+i; e, uj) - Id) - F(-, e, a;)(//+(-, e, u) - Id)), (5.16) 

where I0 is defined by (4.13). In addition, using (5.14), (5.15) formula (5.16) can be 
rewritten as 

= /e(5"F(-, 0, a;)(//+(-, 0, to) - Id) + F"(-, 0, a;)5>+(-, 0, to)). (5.17) 
The following estimates hold: 

||5"F(-,^,a;)(//+(-,^,a;)-/d)||^„,,^^3)< 
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||F-(-,^,a;)5V(-,^,a;)llc«,.(r)< 

gig2e||<^-/IU^,.(R3) (5.19) 
(1 - q2eY 

under the condition that q2S < 1, where qi, q2 are given by (4.43), 

gs = C3(q!, (t)||x|Ic.i+„,.(R3^, (5.20) 

e is the number of (3.29), Ai, C2, C3 are the numbers of Lemmas 3.1, 4.2, 4.3, ^ e 
Proo/ 0/ ("J.i,?;. Using (3.28), (4.39) and the definiions of F, F", we have that 

Using (4.40), (4.42) and Lemma 4.3 we obtain that 



+00 



||x(/.+ (-, 9,00) - /d)||^„,.(R3) < qiq3sjyq2sy-' = (5.22) 

under the condition that q2e < 1. 

Estimate (5.18) follows from (5.21), (5.22) and Lemmas 3.1 and 4.2. 

Proof of (5.19). Using (4.40) for //+ and for ji^^ we have that 

+00 

F^{; e, a;)5>+(., e, u) = F^{; 9, u) ^ D^eS^wji; 9, u), 
S'^Wj{-,9,io) = Wj{-,9,io) - Wj{-,9,u), 



(5.23) 



where w'j is defined as Wj of (4.40), but with F" in place of F. In addition, one can see 
that 

5^w^{;9,u;) = 5-F{.,9,u), 

5'^Wj+^{;9,uj) = 5''F{;9,uj)D_ew]{;9,uj) + F{;9,uj)D_05''wj{;9,uj), 3 en. ' ' 

In addition, using (5.21), (5.24), Lemma 4.3, estimate (4.42) for and Lemma 3.1 we 
have that 

Wwj+i{-,9,uj)\\^^^^^^s^ < 
C3{a,a)\\5''F{;9,oj)\\^,^^^^^^3^qiqi-^e^+ 

C3(a,a)||F(-,^,c^)||^,+„^^^^3^||5"«;,(.,6',a;)||^„_,^^3^ < (5.25) 
qiqie^\6''f\\^^^^^^3^+q2e\\6''wj{-,9,uj)\\^^^^^^3y 
||5-«;i(.,^,a;)||^_^^3^ < gi||5-/||^^,^^^3^. 
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Using (5.25) we obtain that 



Zj+i = {q2ey + q2eZj, Zi = 1. 
In addition, one can see that 

Zj=j{q2ey-\ jeN. (5.27) 
Estimate (5.19) follows from (5.23), (5.26), (5.27), Lemmas 3.1 and 4.3 and the formula 

+ 00 ^ 

= ^, 0<r<l. (5.28) 

Using (5.12), (5.18), (5.19) and Lemma 4.1 we obtain that 
||5-T(-,^,a;)||^o,.(x^)< 

^l^:^+ (l-g2e)^ >'^l'^ /IIl^.^(R^)' ^e^- '^eS- 
Proceeding from (5.29) and assuming that £i of (3.29) is so small that 

92£i < 1/2 (5.30) 

we obtain that 

ll'^"r|lLoc..(A) < (2AiC29i93 + 49i92)e||<^VIUo...(R3^, (5.31) 

where gi, g2, <?3, Ai, C2 are the same that in (5.18), (5.19). 

Finally, (3.31) follows from (5.12), (5.31), (4.47) and Proposition 3.1. 
Proposition 3.3 is proved. 

6. Final remarks 

Remark 6.1. In a subsequent paper we plan to generalize the iterative approach of 
the present work to the case of the polarization tomography with limited phase measure- 
ments, see, for example, [HL], [Sh3] for more information on this problem. Actually, in the 
framework of the model described by (1.1) the polarization tomography with limited phase 
measurements is reduced to the inverse problem for (1.10) with / — {l/2)tr {ngfn0)ld in 
place of /. In this inverse problem we do not plan to restrict the "scattering" matrix S to 
its element only (in contrast with results of the present work). 

Remark 6.2. It remains unclear whether a version of the Riemann-Hilbert prob- 
lem method of [MZ], [N] can be used for solving Problem 1.1, instead of the iterative 
approach of the present work. The reason is that the dependence of F{x,9,(jj) on the 
spectral parameter 6 (and, more precisely, the quadratic dependence on 9 of the element 
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F22{x,6 ,uj) = Q-^ f{x)6-^) in (1.10) is not appropriate, in general, for direct applications of 
the Riemann-Hilbert problem method of [MZ], [N]. 

Remark 6.3. On the other hand (with respect to Remark 1.2), if / is skew-symmetric, 
fij — —fji, then F22 = and the dependence of F{x,9,u) on 9 is appropriate for direct 
applications of aforementioned Riemann-Hilbert problem method to the inverse problem 
for equations (1.1), (1.10). We remind that some results on the polarization tomography 
with skew-symmetric /, including examples of transparent /, were given in [NS]. However, 
the Riemann-Hilbert problem method was not yet used in these studies. 
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